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Abstract. In the framework of the expected association between gamma-ray 
bursts and gravitational waves, we present results of an analysis aimed to search 
for a burst of gravitational waves in coincidence with gamma-ray burst 050915a. 
This was a long duration gamma-ray burst detected by Swift during September 
2005, when the Virgo gravitational wave detector was engaged in a commissioning 
run during which the best sensitivity attained in 2005 was exhibited. This 
offered the opportunity for Virgo's first search for a gravitational wave signal 
in coincidence with a gamma-ray burst. The result of our study is a set of 
strain amplitude upper-limits, based on the loudest event approach, for different 
but quite general types of burst signal waveforms. The best upper-limit strain 
amplitudes we obtain are h rss = O(10" 2() ) Hz" 1 / 2 around ~ 200 - 1500 Hz. 
These upper-limits allow us to evaluate the level up to which Virgo, when 
reaching nominal sensitivity, will be able to constrain the gravitational wave 
output associated with a long burst. Moreover, the analysis here presented plays 
the role of a prototype, crucial in defining a methodology for gamma-ray burst 
triggered searches with Virgo and opening the way for future joint analyses with 
LIGO. 



PACS numbers: 95.55.Ym: Gravitational radiation detectors; 95.85.Sz: 
Gravitational radiation; 97.60.-s: Late stages of stellar evolution; 98.70.Rz: 
gamma-ray sources, gamma-ray bursts 
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1. Introduction 

Gamma-Ray Bursts (GRBs) are intense flashes of 7-ray (and X-ray) photons, lasting 
from few milliseconds to several minutes, followed by a fainter and fading emission 
at longer wavelengths called the "afterglow" [TJ [5] . GRBs are detected at a rate of 
about one per day, from random directions in the sky. They fall into two apparently 
distinct categories, namely short-duration (nominally, less than 2 s), hard-spectrum 
bursts (short GRBs) and long-duration (greater than 2 s), soft-spectrum bursts (long 
GRBs) [HI IH [5J EH ■ Of course, this separation is not strict and the two populations 
do overlap, but such a distinction has suggested two different types of progenitors. 
Progenitors of long GRBs are thought to be massive, low-metallicity stars exploding 
during collapse of their cores; mergers of neutron stars (NSs) probably represent the 
most popular progenitor model of short GRBs at the present time (TJ [2] . 

GRBs are likely associated with a catastrophic energy release in stellar mass 
objects. The sudden emission of a large amount of energy in a compact volume (of the 
order of tens of kilometers cubed), leads to the formation of a relativistic "fireball" 
of pairs, 7-rays and baryons expanding in the form of a jet, while part of the 
gravitational energy liberated in the event is also converted into gravitational waves 
(GWs) [H [2]- In the standard fireball model, the GRB electromagnetic emission is 
thought to be the result of kinetic energy dissipation within the relativistic flow, taking 
place at distances greater than ~ 10 13 cm from the source [UH]. The electromagnetic 
signal can give indirect but important information on the progenitor's nature (e.g. its 
properties can constrain the structure and density of the circumburst medium, and it 
allows the identification of host galaxies). However, to reach a clearer understanding of 
the phenomenon, one should search for a direct signature of the progenitor's identity, 
which may be observed through the gravitational window. The energy that is expected 
to be radiated in GWs during the catastrophic event leading to a GRB would, in fact, 
be produced in the immediate neighborhood of the source. Thus, the observed GW 
signal would carry direct information on the properties of the progenitor. 

In this paper, we present an analysis of the Virgo data [Z] simultaneous with the 
long GRB 050915a [8], with the goal to constrain the amplitude of a possible short 
burst of GWs associated with this GRB. At the time of GRB 050915a, Virgo was 
engaged in a 5-day long data run, named C7. Virgo's sensitivity during C7 exceeded 
that of all its previous runs. The lowest strain noise was ~6x 10~ 22 Hz -1 / 2 around 
~ 300 Hz. This is the first time a study of this kind is being performed on Virgo 
data, so the work presented here aims also to define a procedure of analysis for GRB 
searches with Virgo. In the very near future, these kinds of studies will take advantage 
of the joint collaboration with LIGO, and the existence of an established procedure is 
fundamental. 

The sensitivity of Virgo during C7 was comparable to that of LIGO at the 
time when a coincidence search with GRB 030329 was performed [9], and the 
upper-limits that we set for GRB 050915a are of the same order of magnitude, 
i.e. O Hz~ 1/2 ). The LIGO results on GRB 030329 thus represent a natural 

comparison for our analysis, even if procedures followed in this present study were 
developed for a single-detector search, while those of the LIGO study were for a 
double detector search. 

In what follows, we first present an overview of the Virgo detector and its status 
during C7 (section [2|) . Then, we recall the scenarios for GRB progenitors and their 
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associated GW emission (section [3]) , and we describe the main properties of GRB 
050915a (section^}. Furthermore, we present the analysis of Virgo data in coincidence 
with GRB 050915a (section [3]) , and finally discuss our results (section . 

2. The Virgo detector 

The Virgo gravitational wave detector, jointly funded by INFN (Italy) and CNRS 
(France), is located near Pisa, at the European Gravitational Observatory (EGO). It 
is a power recycled Michelson interferometer with / = 3 km long arms, each containing 
a Fabry-Perot cavity (see e.g. [10] for a recent review of Virgo's status). 

GW interferometric detectors like Virgo [11] have different types of source targets 
for their searches. These can be usefully divided in four main classes: stochastic 
waves, bursts, coalescing binaries and periodic waves. Sources that contribute to the 
first class are e.g. binary stars and primordial GWs (e.g. [T2] [T31 HI]). Sources 
expected to produce GW signals of the other three classes are e.g. compact binaries 
and their coalescence (e.g. [15] and references therein), rotating NSs with a non 
axis-symmetric mass distribution along the rotation axis or NSs instabilities (e.g. 
[El [TTl [18] ) , collapse of massive stars and supernova (SN) explosions (e.g. [19] and 
references therein). GRBs, short and long, are thought to be linked respectively to 
the coalescence of compact binaries and collapse of massive stars [2], and this has 
motivated searches for GWs signals in association with these sources [9, 21-37]. 

By the beginning of September 2005, Virgo was engaged in a commissioning run 
named C7, with the aim to test the gain in sensitivity after several improvements 
were performed on the detector (see [37] for a review of Virgo status during C7) . The 
run lasted 5 days, with a duty cycle of ~ 65%, and a mean sensitivity such that an 
optimally oriented 1.4 - 1.4 M Q NS - NS binary coalescence, at a distance of ~ 1 Mpc, 
would have been detected with a signal-to- noise ratio of 8. The noise spectrum of 
the detector during C7 may be roughly divided into two intervals. Below 300 Hz, 
the main contributions were from the noise associated with the control systems of the 
arms optical axes and of the angular degrees of freedom of the mirrors. Above 300 
Hz, the dominant contribution was from shot noise [38) . Figure [T] shows a comparison 
of the Virgo sensitivity curve during C7, and the LIGO sensitivity during its second 
science run (S2), when a search for a burst of GWs in coincidence with GRB 030329 
was performed [pj. 

The Virgo data acquisition system acquires the interferometer signals, sampled 
at 20 kHz, and a large number of auxiliary signals pertaining to components of the 
interferometer or of its physical environment. This information is used to assess the 
state of various interferometer sub-systems (e.g. the mirror suspension status) and, in 
general, the quality of the collected data. A set of flags summarizes the data quality 
and, in particular, the "Science mode" flag indicates that the interferometer is in a 
stable locked configuration, and the modulation lines used for calibration purposes are 
properly set. The main physical information for the detection of GWs is extracted 
by reconstructing the interferometer strain. This fundamental step consists of the 
extraction of the arms' length relative difference 61/1, i.e. the amplitude of the GW 
signal, from the output dark fringe signal [59] . The last is corrected for the effects 
of controls by subtracting the contribution of the pendulum motion of the suspended 
mirrors, caused by injected control signals on the reference mass coil drivers. Since 
the reconstruction procedure performs a subtraction of the control signals, it removes 
part of the injected control noise and tends to cancel the calibration lines, like any 
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Figure 1. The Virgo sensitivity during C7 run and the Virgo nominal sensitivity 
are plotted in black. Typical LIGO Hanford sensitivities during the S2 arc shown 
in red (2 km) and blue (4 km) |40| . 



signal injected in the mirror control. In practice, the goodness of the calibration 
lines subtraction is used to monitor the quality of the reconstruction. During the C7 
run, the error in the /i-reconstructed data (h(t) = 81/1) was estimated to be around 
~ +20% — 40%. Hereafter, this is assumed as a systematic error in our analysis. 

3. GRB progenitor models and the expected GW signal 

The actual favored scenario for long GRB progenitors is the so-called "collapsar" 
model, which invokes the collapse of a massive star down to a black hole (BH) 
with formation of an accretion disk, in a peculiar type of SN-like explosion (see e.g. 
[HI 1321 S3] ) • On the other hand, the favored scenario for short GRB progenitors is 
the compact binary (NS - NS or BH - NS) merger. This process is believed to occur 
extremely quickly and be completely over within a few seconds, naturally accounting 
for the short nature of these bursts. Unlike with long bursts, there is no conventional 
star to explode and therefore no SN is expected. To produce a GRB, both long or 
short, it is required that the progenitor stellar system ends as a rotating BH and a 
massive disk of matter around it, whose accretion powers the GRB ultra- relativistic 
fireball in the form of a jet, along the rotational axis of the system [JJ. Due to the 
relativistic beaming effect, only observers located within the jet opening angle are 
able to observe the emission from the jet. In the standard assumption (see e.g. [33]), 
the jets are uniform, with sharp cut-offs at the edges, and the line of sight cuts right 
across the jet axis, i.e. the Earth is near the center of the 7-ray beam. 

GWs are expected to be emitted in association with both long and short GRBs 
[351 US]- Coalescing binaries, thought to be associated with short bursts, are one of the 
most promising GW sources for interferometric detectors like Virgo. For such systems, 
a chirp signal should be emitted in GWs during the in-spiral, followed by a burst-type 
signal associated with the merger and subsequently a signal from the ring-down phase 
of the newly formed BH. The last, initially deformed, is expected to radiate GWs 
until reaching a Kerr geometry [46] . In the collapsar scenario, relevant for long GRBs, 
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the high rotation required to form the centrifugally supported disk that powers the 
GRB, should produce GWs emission via bar or fragmentation instabilities that might 
develop in the collapsing core and/or in the disk [46 j. Moreover, asymmetrically in- 
falling matter is expected to perturb the final BH geometry, leading to a ring-down 
phase [46]. While an axis-symmetric rotating collapse and core bounce would give no 
contribution to GW emission along the GRB axis, bar and fragmentation instabilities 
are all dominated by modes with spherical harmonic indices I — m = 2 47J, implying 
that GRB progenitors would emit more strongly along the GRB axis than in the 
equatorial plane (i.e. the orbital plane of the disk fragments or of the bar). The same 
holds for the ringing BH [47]. Thus, in the standard scenario having the Earth near 
the center of the 7-ray beam, the detector is located in the maximum of the emission 
pattern for GWs dominated by spherical harmonic indices I = m = 2. 

Simulations of GW emission in the process of core-collapse in massive stars 
(see e.g. [TO] and references therein), has advanced much recently. However, their 
applicability to GRBs is not clear, since in these simulations the final outcome of the 
collapse are generally NSs. Moreover, those simulations typically do not arrive at the 
point when non-axis-symmetric instabilities are developed. Recently, [48] have further 
extended the post-bounce evolution of some pre-supernova models characterized by a 
relatively quick spinning stellar iron core, finding that these models do indeed develop 
a considerable non-axis-symmetry after the bounce, with associated GW emission in 
agreement with the expectations from a spinning-bar model [48]. The study of GW 
emission by gravitational collapse of uniformly rotating NSs to BHs, is also advancing 
(see e.g. [49]). However, due to the high degree of axis-symmetry of these simulations, 
their outcomes do not represent the best scenario for GRB progenitors. In this article, 
we assume a model based on a best case scenario, and work under the hypothesis that 
bar and fragmentation instabilities do indeed develop, and that the Earth (i.e. the 
Virgo detector) is aligned with the center of GRB 050915a 7-ray beam. 

A final important aspect is the expected delay between the electromagnetic trigger 
and the associated GW signal. Typically, a GW signal is searched for within a window 
of 180 s around the GRB trigger. For a GW burst associated with the formation of 
the GRB central engine, i.e. with the launch of the jet that subsequently powers the 
GRB, the electromagnetic trigger should follow the GW trigger. As described in [50] . 
in the case of a collapsar (relevant for our analysis), the time delay between the two 
triggers is dominated by the time necessary for the fireball to push through the stellar 
envelop of the progenitor, which can be of the order of 10 — 100 s. Thus, a period 
of 120 s before the trigger time is typically selected to over-cover these predictions. 
Moreover, given that some models predict a GW signal contemporaneous with the 
GRB emission (i.e. extending from the time of onset of the GRB to its end, see [5"T]). 
the 60 s after the electromagnetic trigger are also included in our search. 

4. GRB 050915a 

On the 15 th of September 2005, at T= 11 : 22 : 42 UT, the "Burst Alert Telescope" 
(BAT) on-board Swift [52] triggered and located GRB 050915a [SJ. The BAT on- 
board calculated position was RA=05h 26m 51s, Dec=-28d 01' 48" (J2000), with 
an uncertainty of 3 arcmin. The BAT measured a T90 [jj duration of 53 ± 3 s in the 

I The Tgo duration is defined as the time necessary to collect from 5% to 95% of the total counts in 
the specified energy band. 
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15 — 350 keV energy band [53], thus GRB 050915a was classified as a long- type GRB. 
The "X-Ray Telescope" (XRT) began observing the BAT position at 11 : 24 : 09 UT 
(~ 87 s after the trigger, [S3]). The new refined position was RA=05h 26m 44.6s, 
DEC=-28d 01m 01.0s [53J. Finally, the Swift "Ultra- Violet and Optical Telescope" 
(UVOT) began observing the field of GRB050915 - 85 s after the BAT trigger [f>5], 
and the optical and IR follow-up of this burst was performed by different telescopes 
[561 EH [58] (SSI HE EU- In the radio band, VLA observations on September 18.58 UT 
revealed no radio source in the error circle [62] ■ Recently, evidence has been found for 
a possible distant and/or faint galaxy [53] [M], however the redshift of this burst still 
remains unknown. 

5. Coincidence analysis for a burst of GWs 

Given the lack of accurate predictions on the expected GW waveforms that might 
be associated with long GRB progenitors, we have chosen in this search to look for 
a GW burst- type signal associated with GRB 050915a in a model independent way. 
Our analysis aims to: 

• set-up a procedure for GRB triggered searches with Virgo, that goes well beyond 
the single case of GRB 050915a, but has a much more general interest, both in 
the context of GWs triggered searches and GRB research, and also in view of the 
future joint collaboration with LIGO; 

• constrain the amplitude of the associated GW emission to quantify Virgo's 
capability to provide information pertaining to models for GW production by 
GRBs, and also to define possible margins of improvement in view of the expected 
enhancement in sensitivity; 

• define an algorithm that will find application also for the categories of short GRBs 
(merger and ring-down phase, see section[3J, which are expected to be nearer and 
thus more promising sources of GWs than long bursts. 

5.1. Wavelet Analysis 

For our analysis, we relied on a new wavelet-based transient detection tool, the Wavelet 
Detection Filter (WDF). Wavelets were introduced in the 80's as a mathematical tool 
to represent data both in time and frequency [65]. The wavelet transform is defined 
as the correlation of the data x(t) against the wavelets tp a ,b, 



The wavelet family is obtained by translations and dilations of a reference waveform 



The wavelet transform gives a representation of the signal in terms of the scale a 
(associated to frequency) and time b. The reference wavelet ip(t') is chosen to be a zero- 
mean function of unit energy, well-localized both in time (around t' = 0) and frequency. 
Consequently, this is a short-duration waveform with few cycles. Wavelet-based 
representations are well-suited for burst-like signals because of similarities between 
those signals and the analyzing wavelet ip. 




(1) 




(2) 
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The wavelet family is redundant but when the sampling of the scale a and time 
b axes is dyadic i.e., when (aj,bk) = (2 J , fe2 J ) for j > and k integers, it forms an 
orthonormal basis, provided some geometrical constraints on the choice of ip are set 
(we will not detail them here, the reader is referred to [55]). The wavelet transform 
produced in this way is referred to as a discrete wavelet transform, and we use it to 
analyze the data. Since we are dealing with signals sampled at a rate of f s = 20 kHz, 
and we consider data blocks of N samples, the dyadic sampling is limited to the range 
j = 0, . . . , log 2 N and b k /aj = k/f s with k = 0, . . . , N - 1. 

5.2. Best matching wavelets and thresholding 

The wavelet transform can be viewed as a bank of matched filters. We select the best 
matching wavelets (largest correlation coefficients) by thresholding. Let us define the 
soft-thresholding operator T(w; r)) = sign(w)(\w\ — rf) if |to| > 77 and otherwise, and 
the thresholded coefficients of the discrete wavelet transform Wj t k = T(W x (aj,bi c );rj). 
We use the Signal to Noise Ratio (S w ) as a statistic 



to discriminate the presence or absence of a burst-like signal in the data, with being 
an estimate of the variance of the noise. The threshold choice is 77 = y/2 log Na n . 
This choice is linked to a general result by Donoho & Johnstone |66j concerning non- 
parametric denoising. Let our data x(t) = s(t) + n(t) be the sum of a signal s(t) and 
Gaussian white noise n(t). It can be shown that the signal estimate obtained as 



i.e. by inverting the thresholded wavelet transform, minimizes the mean-square error 
over a broad class of signals [67] . 

5.3. Pipeline: preprocessing, analysis and trigger selection 

We preprocess the data by applying a time-domain whitening procedure [BP] using an 
estimate of the noise spectral density given by an autoregressive (AR) fit over the first 
1000 s of the data set. We divide the whitened time-series into overlapping blocks 
of duration Atw — 12.8 ms, which also sets the time resolution of the search. The 
epochs of two successive blocks differ by At s = 0.6 ms. 

We already mentioned that the wavelet basis may be interpreted as a template 
grid. It is well-known that redundant grids are better suited for burst detection 
because they increase the chance of a good match between the signal and one of the 
templates. Instead, the wavelet bases are sparse by construction. To compensate 
for this sparsity, we insert redundancy by combining the results from several wavelet 
decompositions. We empirically select 18 different wavelets (including Daubechies 
wavelets from 4 to 20 [BS], the Haar wavelet and the windowed Discrete Cosine [67 ). 
For each data block, we compute these 18 wavelet transforms. Thus, we obtain 18 
S w estimates and select the largest one (SNR). Furthermore, we generate a shorter 
trigger list by selecting only SNR values larger than 4. 

Typically, the energy content of a burst signal will be tracked by the WDF in 
the form of a "cluster" of successive triggers, each containing a different fraction of 




(3) 




(4) 
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Figure 2. BAT (15-350 keV) light curve of GRB 050915a in total count rate 
(counts/s, thin solid line). Data of the BAT light curve for GRB 050915a has 
been downloaded from the online archive 1681 . The thick solid line defines the 
time length and position of the on-source region, the thin-dashed segments mark 
the 60 s of data before the start and after the end of the signal region which 
are excluded from the analysis (so to separate the background from the signal); 
the thick-dashed segments mark the portions of the off-source region around the 
GRB trigger time. The whole off-source region considered in our analysis is much 
longer, extending on the left of the plot up to —9973 s, and on the right of the 
plot up to 6863 s. 



its energy. To assign to each candidate event a unique time and SNR value, such 
clustering should be properly organized. As such, we cluster all triggers having a time 
difference less than 10 ms (~ Atw)- For each cluster, a candidate event is defined 
having the time and SNR of the trigger at which the maximum SNR of the cluster 
is reached, i.e. where the signal leaves most of its energy. 

5.4- The analysis method 

For our analysis, we rely on a single stretch of data during which the interferometer 
was maintained in the same configuration and the /i-reconstruction processes was 
recognized as good. Such a stretch is between the GPS times 810808602 s and 
810825638 s, for a total of 16836 s, containing the GRB trigger time (i.e. GPS time 
~ 810818575 s). We define as the on-source region a data segment 180 s long, 120 s 
before the GRB trigger time and 60 s after (see the thick-solid line in Fig. [2|). This is 
the time window where we searched for a coincidence with the GRB trigger. The rest 
of the data in the stretch, with the exception 60 s before the start and after the end 
of the signal region (see the thin-dashed in Fig. [2]), belong to what we define as the 
off-source region (see the thick-dashed lines in Fig. [2]). The analysis of the off-source 
region is used to assess the data quality and to study the statistical properties of the 
background. As explained in section [31 the on-source region has been chosen to start 
120 s before the trigger time, so as to over-cover most of astrophysical predictions 
regarding the expected delay between the GRB and the associated burst-type GW 
signal. Moreover, in view of models that predict a GW signal contemporaneous with 
the GRB emission and considering that GRB 050915a had a Tgo duration of 53 ± 3 s 
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(see section 2]), we have chosen our signal region to end 60 s after the trigger. The 
same choice for the time-length of the on-source region was also implemented in [S], 
for the case of the long GRB 030329. 

Our pipeline is calibrated by adding simulated signals of various amplitudes and 
waveforms to data in the off-source region. The simulated signals were produced using 
the Virgo SIESTA simulation code [70] . The resulting data stream is processed in the 
same way as for the off- and on-source regions. Using simulated signals we evaluate the 
detection efficiency as a function of the simulated signal strength, which we quantify 
in terms of root-sum square amplitude of the incoherent sum of the contributions from 
the "plus" and "cross" polarizations: 




h%(t))dt. (5) 



This allows us to make physical interpretations with the results. 

The calibration procedure is based on simulations of plausible but quite general 
burst-type waveforms, with different amplitudes, characteristic frequencies and 
durations, chosen on the basis of the considerations explained in section 15.51 The 
times at which those signals are added to the off-source data are randomly determined 
by following a Poisson distribution with a mean rate of 0.1 Hz. 

Knowledge of the source position is also used when adding the simulated signals 
to off-source data, by considering the antenna response at the GRB position and 
time. A GW arriving at the interferometer from the GRB direction can be described 
as a superposition of two polarizations amplitudes h + and h x . The response of the 
interferometer to such a wave is given by [71) : 

Sl/l = h{t)=F+h+(jt) + F x h x (t) (6) 

where F + and F x are expressed as functions of the source position and of the wave 
polarization angle tp |72j . The antenna patter functions F+ and F x can be written as: 

F + = F° cos(2?/>) - F x sin(2V>) (7) 

F x = Fl sin(2^) + F° cos(2^) (8) 

where F + (ip = 0) = F° and F x (tp = 0) = F°. In the case of GRB 050915a, we have 
F° ~ 0.32 and F° ~ 0.21. 

5. 5. Choice of plausible waveforms 

5.5.1. Gaussian waveforms To calibrate our pipeline, considering the great 
uncertainties in the waveforms associated with long GRB progenitors, we added 
different burst-type signals to the off-source data. Our simplest choice was for 
Gaussian signals, having the following form: 

h(t) = h exp[-(t - t ) 2 /2a 2 ]Fl (9) 

where to is the time at which the signal is added to the off-source data stream, and a 
values of 0.5 ms, 1 ms and 1.5 ms were considered. These were broad-band, linearly 
polarized waveforms along the + direction, with the unknown polarization angle ip 
set to zero. Note that for a given value of a and ho, a different choice of ip rescales 
the waveform amplitude arriving at the detector by a factor of cos(2t/>), while leaving 
unchanged its shape over the detector bandwidth. Thus, the efficiency curves for the 
general ip case can be estimated from the ip = ones here presented, by rescaling the 
h rss corresponding to a given detection efficiency for a factor of COB L^,s . 
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5.5.2. Sine- Gaussian waveforms To mimic GW emission by GRB progenitors during 
the phase of collapse, fragmentation or bar instabilities, we considered sine-Gaussian 
waveforms. Taking a best case model scenario of GW emission from a triaxial ellipsoid 
rotating about the same axis as the GRB (i.e., the direction to the Earth, see Eqs. 
(|A.8P and (|A.9[) in the appendix | Appendix A[ ), and using a Gaussian amplitude as 
the simplest way to mimic the impulsive character of a GW burst, we consider signals 
having in the wave frame the following form: 

h+ = h exp[-(t - t ) 2 /2a 2 ] cos(2tt f {t - i )) (10) 
h x = h Q exp[-(t - h) 2 /2a 2 } sin(27r/ (t - t Q )) (11) 
for an unknown value of polarization angle The detector response to such types of 
signals is then computed using Eqs. CO)-©. After some algebra, the resulting h(t) 
can be written as: 

h(t) = h exp[-(t - to) 2 /2a 2 ][Fl cos{2tt f Q (t - h) - 2ip) + 

+F°sin(27r/o(t-io)-2^)] (12) 

In our analysis we set if) = and we span the frequency range fo ~ 200 — 1500 Hz, 
as suggested by the predictions for GW emission from GRB progenitors, when 
fragmentation or bar instabilities arc developed (see dash-dotted lines in Figs. 3-5 
of [46]). For each fo, we consider two values of Q, i.e. Q — 5 and Q — 15. 

It is worth noting that, for signals of the form (110p - ([Ti"]) with Q — 2irfoa > 3 (i.e. 
for relatively narrow-band signals), one has: 




(13) 



and: 



h?{t)dt ~ h rss \ 1 + > . 1 x> (14) 



where h(t) is given by Eq. (|12|). I n this approximation, if the detector noise is 
roughly constant within the relatively narrow signal bandwidth, the detected SNR is 
proportional to the above integral. Thus, the detection efficiency as a function of the 
h rss is expected to be independent on the choice of ip. 



5.5.3. Damped sinusoid waveforms Consider now the phase of BH ringing. A Kerr 
BH distortion can be decomposed into spheroidal modes with spherical-harmonic-likc 
indices I and m (see e.g. [73]). The quadrupole modes (I = 2) presumably dominate 
[73], while the paramount m- value depends upon the matter flow. In particular, the 
m = ±2 modes are bar-like, co-rotating (m = +2) and counter-rotating (m = —2) 
with the BH spin, and the I = m = 2 mode is expected to be the most slowly 
damped one. As underlined in [74J, numerical simulations of a variety of dynamical 
processes involving BHs show that, at intermediate times, the response of a BH is 
indeed well described by a linear superposition of damped exponentials. Generally 
speaking, the polarization of the ring-down waveform will depend on the physical 
process generating the distortion of the BH (see e.g. [75]). Since the I = m = 2 mode 
may be preferentially excited in the presence of binary masses or fragmentation of 
a massive disk, it is commonly assumed that the distribution of the strain between 
polarizations h+, h x for this mode mimics that of the in-spiral phase [47l [76l [74] . 

h+ = ho-{l + cos 2 d ) cxp(-t/r) cos(2vr/ i + £)8(t), (15) 
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h x = /i cos 6» exp(-t/r) sin(27r/ t + £)©(*), (16) 

where £ is an arbitrary phase, 9q is the inclination of the angular momentum axis with 
respect to the source direction in the sky, and Q(t) is the normalized step function. 
Since we expect to observe the GRB on-axis, this polarization is also circular. 
Thus, we added to the off-source region signals having the form 

h(t) = h exp[-(t - i )/r] x 

x [F° cos(2^/ (t - t ))Q (l - —L—] 0(t - t ) + 



4/o(*-*o). 

+F°sm(2irf (t-t ))Q(t-t )], (17) 

where again 9(t — to) and 6 ^1 — 4^-^) ^ are normalized step functions 5J 

The characteristic frequency of the I = m = 2 quasi-normal mode of a Kerr BH 



is estimated as 



f = 32 kHz (1 - 0.63(1 - a) 3 / 10 ) (J£j , (18) 

where a is the dimensionless spin parameter of the Kerr BH, while the damping time 
can be estimated as 

r^(AJ)- 1 =nf /Q(a), (19) 

where Q(a) = 2(1 - a)" 9 / 20 @g [77]. For GRB progenitors it is typically assumed 
a = 0.98, since the BH is supposed to have spun up to near maximal rotation by a 
massive accretion disk 46 . In the collapsar model, M is expected to be of the order 
of ~ 1 Mq, which implies /q ~ 10 kHz [46], so a detection would be difficult. 

However, the collapsar model for long GRBs is one of a larger class of proposed 
progenitor models, all leading to a final BH plus accretion disk system. Thus, the 
process of GW emission can always be described in a way similar to the collapsar 
case |46j : a high rotation rate producing bar or fragmentation instabilities in the disk, 
followed by a BH initially deformed encountering a ring-down phase. Among these 
variants of the collapsar model, the BH-white dwarf scenario may be characterized by 
higher BH masses (M ~ 10 Mq) and lower /q values, down to ~ 1 kHz (see Fig. 3 
of [46) ) . According to Eq. (|19p . a fo around 1 kHz would imply a damping time of 
~ 0.3 ms. Thus, we choose to span a frequency range between ~ 800 Hz and 3 kHz, 
for r values of 0.3 ms, 1 ms and 1.5 ms. 



6. Results and Discussion 



6.1. Data Quality 

We applied the WDF to the data of the off-source region and we derived the 
distribution of the trigger strength. In the off-source region, we selected ~ 2.1 x 10 3 
triggers crossing the threshold SNR = 4. We processed the resulting list so as to 
eliminate triggers related to instrumental artifacts. The trigger rejection operates in 
two steps [78], namely a data preselection followed by a glitch removal procedure, that 
we can summarize as follows. 

K The reason for multiplying cos(27r/o(t — to)) by ^1 — 4f (t—t ) ) ' s *° av °id a discontinuity at the 
beginning of the waveform, which would result in an infinite energy, even though h rss would remain 
finite. 
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Figure 3. The mean SNR distribution found in the off-source region (black) and 
its ±lrr interval (blue) is compared with the on-source SNR distribution (red). 
As is evident, the on-source distribution is within the lrr interval around the 
mean off-source one, confirming that the on-source and off-source distributions 
are statistically compatible. 



First, triggers falling into periods where known instrumental problems occurred 
(e.g. saturation of the control loop electronics, problems in the /i-reconstruction 
process) or when aircraft (known to produce high seismic/acoustic noise which couples 
into the interferometer) fly over the instrument, are discarded. This preselection, while 
leaving the on-source region untouched, cuts out from the off-source segment (16500 s) 
about 14 s associated with high acoustic noise, occurring about 1945 s after the start 
of the off-source stretch. Furthermore, a glitch removal procedure was applied. An 
extensive study has been performed to establish the correlation of triggers produced 
by the burst search pipelines, and environmental or instrumental glitches occurring 
during C7 (see [78] for a detailed description). This study provided the definition of 
a series of veto criteria, based on information given by the auxiliary channels, and 
introduced a "dead time" of ~ 6.3% [78 on the complete data set of the C7 run. The 
application of these criteria in our analysis leads to a dead time of ~ 582 s in the 
off-source region (i.e. ~ 3.5% of its duration), and of ~ 1.3 s in the on-source region 
(i.e. ~ 0.75% of its duration). The loudest on-source event remains unaffected by the 
veto procedure. 

It is important to note that we have not applied all the vetoes designed for 
the C7 data. In fact, the majority of spurious burst triggers has been related to 
the so-called "Burst of Burst" (BoB) [75]. BoBs originate from a misalignment of 
the interferometer mirrors, which increases the coupling of the laser frequency noise 
into the interferometer itself and causes a noise increase lasting up to a few seconds. 
Procedures have been defined to veto the BoBs. However, as estimated from the 
complete dataset of C7 run, these introduce a large dead time (~ 16% of the run 
duration [75])- Thus, since the BoB cuts can significantly affect the integrity of the 
on-source segment, we chose not to apply them. It is worth noting that during a BoB, 
the interferometer is sensitive to GWs, therefore we can still observe a possible GW 
counterpart of a GRB, albeit with a lower SNR. 
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Figure 4. Distribution of event strengths in the on-source region (number of 
events vs detected SNR). 



6.2. Statistical Analysis 

We show in Fig. [3] the SNR distribution found in the off-source and on-source regions, 
where we have applied the data quality cuts described in section 16.11 The SNR 
distribution in the on-source region is confined below SNR = 9 (see Fig. @|. To test if 
the distribution of events observed on-source is compatible with being only noise, we 
performed different checks. First, starting from the beginning of our data stretch (i.e. 
GPS 810808602 s), we sampled our background distribution by dividing the off-source 
region in ~ 90 successive windows 180 s long. We find that the percentage of such 
windows having a loudest event with SNR > 9 is of about 89% (after applying in each 
of the windows the cuts described in the prevoious section). Second, we computed the 
mean SNR distribution on the 90 off-source windows, and derived for each SNR bin 
the corresponding a. As shown in Fig. [3l the on-source distribution is well within the 
±lcr interval around the mean off-source one, thus being compatible with noise. 

From these tests we conclude that the on-source events are consistent with noise 
and that no clear evidence is found for an exceptional event, with respect to the 
background statistics, that could possibly be associated with the GRB. Thus, we 
move to the definition of an upper- limit, by following the procedure described in 
|79j . To this end, we use simulated signals to determine the strain necessary to have 
90% frequentist probability for such signals showing up as events with SNR > 9, 
i.e. with a SNR above that of the loudest event observed in the on-source region. 
This means estimating the efficiency e at which the instrument and filtering process 
can detect burst events with SNR > 9. According to [79j, the simulated signals are 
added to the off-source data, so as to evaluate e with good statistics, thanks to the 
long duration of the chosen off-source stream. Before proceeding with the efficiency 
estimates, basic sanity checks aimed to guarantee the consistency of our approach 
were applied. These included a comparison of the off-source and on-source regions 
in the time-frequency domain, and a Kolmogorov-Smirnov test (90% confidence level, 
two-sided test) between the on- and off-source SNR distributions. 
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6.3. Detection efficiency and upper-limit strain 

The efficiency e in detecting signals with SNR > 9, is estimated by computing for 
each kind of chosen waveform (see section [5~5)l , the percentage of simulated signals 
found by the pipeline, as a function of the injected strain amplitude h rss . A simulated 
signal added to the noise at a given time to, is recognized by the pipeline at tdet if 
\tdet — to\ < 20 ms. This ±20 ms coincidence window takes into account the duration 
of the wavelet decomposition window (12.8 ms), allowing a partial overlap. Moreover, 
a coincidence window of ±20 ms contains the ~ ±2er portion of the longest duration 
simulated signal (sine-Gaussian waveforms with Q = 15 and frequency 203 Hz, having 
a ~ 12 ms). 

In Fig. [5] and Fig. [5] we show the results obtained for sine-Gaussian waveforms 
with Q = 5 and Q = 15 respectively, at different frequencies. As one can see, 
the detection efficiency depends on the signal frequency. There are two main 
elements determining such dependence: (i) the detector noise level and (ii) the signal 
duration with respect to the window in which the wavelet decomposition is performed. 
Concerning point (i), consider two sine-Gaussian signals at different frequencies, with 
equal strain amplitude h rss . Those will be detected at different SNR values, since 
the detector noise level changes with frequency. Thus, the lowest is the detector noise 
around the signal characteristic frequency, the highest will be its detection efficiency. 
Concerning point (ii) , sine-Gaussian waveforms with the same Q but differing /q have 
different durations (i.e. 5er — 5g3g). The filter capability in detecting a signal is 
maximized when the duration of the window in which the wavelet decomposition 
is performed, is comparable to the signal duration. When taking a shorter wavelet 
window, part of the signal power is lost. On the other hand, choosing a wavelet 
window much longer than the signal duration, implies that the wavelet decomposition 
is dominated by the background, resulting in a loss of efficiency. This is the reason 
why in our analysis, as a trade off, we set a wavelet decomposition window of 12.8 ms, 
comparable to the a ~ 12 ms of the longest injected event. 

In Figs. [7][5]we report the results obtained for damped-sinusoid waveforms with 
t values of 0.3 ms, 1 ms and 1.5 ms. For the given h rss and /o values, the detection 
efficiency decreases with decreasing r. On the other hand, for a given r but different 
characteristic frequencies, the detection efficiency decreases with increasing /q. 

Finally, in Fig. [10] we show the results obtained for the simplest type of simulated 
waveforms, i.e. the Gaussian ones. The 5cr duration of these signals are between 1.5 ms 
and 7.5 ms. For a given h rss value, a higher a in the time domain implies that the 
energy of the Gaussian is in the low frequency region of the detection bandwidth, 
hence the detected SNR (and thus the detection efficiency) is lower. This causes the 
detection efficiency to decrease for increasing a. 

To estimate the h rss at 90% pipeline efficiency for signals with SNR > 9, we fit 
the points in Figs. [511101 with sigmoid functions of the form: 



1 + cxp[-(log/i rss -pi)/p 2 ] 

The 90% h rss values reported in tables [T][3] (third column) are obtained by setting 
in Eq. (|20j) p\ and pi equal to their best fit values, and determining the h rss 
corresponding to the 90% efficiency, e*. The reported errors correspond to the ±2<r 
uncertainty on the best fit curve. We stress that the derived upper-limits (last 
column in tables [1|[3]) , are affected by a +20% — 40% systematic error, related to the 
uncertainties in the calibration of h— reconstruction (see also section [2]). The lowest 
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Figure 5. Efficiency in detecting sine-Gaussian signals with Q = 5, for different 
characteristic frequencies /o, as a function of the injected h rss . Only events found 
with SNR > 9 and within ±20 ms the injection time are counted in the measured 
efficiency. The dashed line marks the 90% efficiency level. For each waveform we 
plot the measured efficiency (points), the best fit sigmoid function and the other 
two curves which account for fitting errors (in most cases, these two curves appear 
overlapped on the best fit one). 



Sine Gaussian waveforms 0=15 




TO" 20 10"' 9 

hrss (Hz"") 



Figure 6. Similar to Fig. [5] efficiency curves for sine-Gaussian waveforms with 
Q = 15 (see the caption of Fig. [5] for further details). 
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Figure 7. Efficiency curves for damped-sinusoid waveforms with damping time 
t = 0.3 ms, for different characteristic frequencies /o, as a function of the injected 
hrss (see the caption of Fig. [5] for further details). 




Figure 8. Similar to Fig. [7] efficiency curves for damped-sinusoid waveforms 
with r = 1 ms. 



h rss upper-limit is obtained for the sine-Gaussian waveform at frequency /o = 497 Hz 
with = 5, for which h?g ~ 2.09 x 1CT 20 Hz" 1 / 2 . 

We compare our results for sine-Gaussian waveforms with Q=5, with those 
obtained by :9] for sine-Gaussian waveforms with Q — 4.5, in association with 
GRB 030329 during LIGO S2. The sensitivity of the Hanford detectors during 
S2 was similar to Virgo during C7. Moreover, the visibility of GRB 030329 from 

LIGO (y^) 2 + (T°) 2 = 0.37), was nearly equal to the one of GRB 050915a from 

Virgo (0.38). LIGO's lowest strain upper-limit, h rss = 2.1 x 10 -20 Hz" 1 / 2 (note 
that according to the different definitions, the upper-limits reported in Table I of [9] 

/(F°) 2 + (F ) 2 

should be divided for y — ± — ^ — ~ — before comparing with the results reported in our 
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Figure 9. Similar to Figs. 17181 efficiency curves for damped-sinusoid waveforms 
with t = 1.5 ms. 




Figure 10. Efficiency curves for Gaussian waveforms with different a, as a 
function of the injected h rss (see the caption of Fig. [5]for further details). 



Table [T] where the quoted h rss strains do not contain the attenuation for the antenna 
pattern (see Eq. [5])), was obtained at fo ~ 250 Hz. In our case, we get the lowest 
value of h rss ~ 2.09 x 10~ 20 Hz -1 / 2 at ~ 500 Hz. At higher frequencies, around 
1000 Hz, the LIGO upper-limit is h rss — 6.5 x 10~ 20 Hz^ 1 / 2 , to be compared with 
h rss ~ 2.96 x 10" 20 Hz" 1 / 2 in the Vir go case. We stress that the LIGO procedure is 
based on the cross-correlation between the output of the two Hanford detectors, while 
our search is a single detector analysis. 

6.4- Astrophysical interpretation 

As described in section 15.51 GWs could give a direct information on the GRB 
progenitor's identity. Of course, the critical aspect in theoretical models for the 
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Table 1. h rss upper-limits for sine-Gaussian waveforms (see Figs. I ■">!('[ ). The 
first two columns give details on the waveform parameter state, the third column 
is the h rss for which 90% efficiency is reached in detecting simultaed signals 
at SNR > 9. The error-bars account for the errors on the best fit values of 
the sigmoid function parameters. Note that these h rS s values are affected by a 
systematic error, as described in Sec. [5] and Sec. 1(131 
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Table 2. hrss upper-limits for damped-sinusoid waveforms (see Figs. 17180 . The 
first two columns give details on the waveform. The last column gives the h rBB 
for 90% efficiency. Note that these h rss values are affected by a systematic error, 
as described in Sec. [2] and Sec. 16.31 
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Table 3. h rss upper-limits for Gaussian waveforms (see Fig. HOP . The first 
column gives details on the waveform. The last column gives the h rss for 90% 
efficiency. Note that these h r3s values are affected by a systematic error, as 
described in Sec. [2] and Sec. 16.31 
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production of GWs in association with long GRBs, is the fraction of energy expected 
to be emitted in GWs, Eqw, during the phases when dynamical instabilities develop. 

Sources radiating energy Eqw? could produce an extremely small h(t) signal at the 
detector, depending on the emission pattern. Nevertheless, we can always associate 
a strain h(t) at the detector with some minimum amount of Eqw radiated by the 
source, selecting an optimistic emission pattern. This is in fact the spirit of the 
analysis presented here, where attention was mostly devoted to that phases of GW 
emission dominated by a I = m = 2 emission pattern (i.e. having a maximum along 
the line of sight). Considering that the redshift of GRB 050915a is not known, in what 
follows we will assume that this burst was at the distance of the nearest long GRB ever 
observed, i.e. GRB 980425 at d^ ~ 40 Mpc. Moreover, we will focus attention on the 
h rss upper-limit obtained for the sine-Gaussian waveform at frequency fo = 203 Hz 
with Q — 5. As shown in | Appendix B| in this case the radiated energy is computed 



When Virgo is running at its nominal sensitivity, the noise strain around ~ 200 Hz is 
expected to be about a factor of 15 lower than during C7 (see Fig. []}. Thus, if we 
assume to have a SNR distribution confined below SNR — 9, then the energy upper- 
limit given in Eq. (|21[) would be lowered of a factor of ~ 225. Further improvement 
may also come in the case of optimal orientation: e.g. if GRB 050915a was optimally 
oriented with respect to the Virgo antenna pattern, the upper-limit in Eq. (|22|) would 
be a factor of ((F£) 2 + (F®) 2 ) — 7 lower. Moreover, the joint collaboration with 
LIGO will help in setting upper-limits, since a coincidence search using three or four 
detectors will be a powerful tool in reducing the tail observed in the SNR distribution 
of the on-source region. 

Some of the most optimistic predictions for the emission of GWs when instabilities 
develop in the rotating core of the massive GRB progenitor or in the disk surrounding 
the final BH, give an upper-limit estimate of the order of ~ O.IMq (e.g. [46] for 
the case of a merger of two blobs of 1 Mq each, formed in the fragmentation of a 
collapsing core). We thus conclude that, under the optimistic assumptions of optimal 
orientation and distance of 40 Mpc, the Virgo detector at its nominal sensitivity will 
start reaching the level of theoretical upper-limit estimates for GW emission by long 
GRB progenitors. 



as: 




(21) 



(22) 
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7. Conclusion 

We have presented the first analysis of Virgo data in coincidence with a GRB trigger, 
aimed to search for a burst of GWs associated with the long GRB 050915a, occurred 
during Virgo C7 run. We have analyzed a time window of 180 s around the GRB trigger 
time, and about 4.6 hours of off-source data, corresponding to a single lock stretch. 
The result of this analysis is a set of loudest event upper-limits on the strain of an 
astrophysical GW signal occurring in association with GRB 050915a. The evaluation 
of the pipeline and detector efficiency for detecting signals showing up as events with 
SNR above the loudest observed in the on-source region, was performed by adding 
a set of simulated burst-type signals to the off-source data, at randomly selected 
times. The waveforms of the simulated signals were chosen taking into account present 
uncertainties in the predictions for GW emission associated with GRB progenitors, 
linked with the ringing of the final BH. 

The best upper-limit strain amplitudes obtained in our analysis are of the order of 
h rss = (2 - 4) X 10~ 20 Hz" 1 / 2 around - 200 - 1500 Hz, affected by a - +20% - 40% 
systematic error. On the basis of these results we conclude that, when running at 
nominal sensitivity, Virgo will start putting interesting astrophysical constraints for 
GW emission in association with GRBs at distances comparable to GRB 980425. 

Short GRBs, probably associated with the merger of compact binaries and 
occurring at lower redshifts [801 181j . will represent even more promising targets. The 
procedure for the analysis presented here may in fact also be extended to the study 
of these sources, especially for the phases of merger and ring-down, but also for the 
last stages of the earlier in-spiral phase. The kind of search implemented here, while 
expected to be less efficient than a matched filtering approach, has the great advantage 
of avoiding strong dependence on exact knowledge of the in-spiral waveforms. Finally, 
these kinds of studies will be of great benefit for the joint collaboration with LIGO, in 
view of which is the hope that a coincident detection could occur with three or four 
interferometers, during the explosion of a relatively near GRB. 
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Appendix A. Details on the sine-Gaussian waveform choice 

In what follows, we address the question on how to mimic GW emission from collapse, 
fragmentation or bar instabilities. We review some results reported in the literature 
that are useful for our GRB analysis. 

Consider a source of gravitational radiation characterized by a mass quadrupole 
tensor Dij [82] • The transverse-traceless components of the metric perturbation 
are related to the transverse-traceless components of the quadrupole tensor by the 
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following relation: 

hJJ = ^-/Djjiae (A.i) 

where i,j, k = 1,2,3. The waveforms of the radiation received by an observer at 
distance d (much greater than the dimensions of the source) can be taken to be 
plane and depend on the relative (angular) orientation between the observer and the 
source. In general, using a system of orthonormal spherical coordinates (r, 0, 4>) and 
considering an observer direction making an angle Oq with the cc3-axis and 4>q with 
the xi-axis, the two polarization states in the plane perpendicular to the direction 
of propagation can be characterized by the following expressions for the two non- 
vanishing components of the perturbation to the galilean metric |82j : 

2G 

G 

h+ = he a e a = -/i0 o o = ~^^( D s 6o ~ D 4>o4><>) ( A - 3 ) 

where Z?0 O O , Dg 0< f, , D^, g are the physical components of Dij projected along the 
directions of the spherical unit vectors. 

Consider now the particular case of a system characterized by a mass quadrupole 
tensor: 



D n 


D 12 


















D 33 



D a , = D 21 D 22 (A.4) 
V D 33 J 

with respect to a set of fixed inertial axes (xi, x 2 , x 3 ), where a^-direction could assume 
the invariant one of the angular momentum or the rotation. No such physical meaning 
is assigned to the x\- and £2 -axes. A large class of realistic astrophysical systems turn 
out to have a mass quadrupole tensor of the form given in Eq. (|A.4[) . These include, for 
example, binary systems, rotating ellipsoidal objects and pulsating/rotating ellipsoids 
[53"] . The latter case, for large amplitude pulsation, corresponds to explosion and 
collapse [SH [5S] . The mass quadrupole tensor in Eq. (|A.4[) is thus relevant for the 
case of a GRB progenitor, for which we expect the £3-axis to be the rotational axis 
of the collapsing core (and of the subsequent BH plus accretion disk system), along 
which the GRB jet is launched. 

For such systems, using the freedom of rotation about the £3-axis, so that a 
reference system can always be chosen to have the observer on the X\-x 3 plane (</f>o = 0), 
one has: 

2G 

ft x = ~TT3^12 cos #0 (A. 5) 

3c 4 a 



G 
3cM 



(D n - D22) + (D 33 -Dn) sin 2 9 (A.6) 



Using the simplified assumption of a rigid, uniform, ellipsoid rotating with an 
angular velocity lo around the X3-axis, the quadrupole mass tensor has a time- 
independent expression in the frame x[ ,x' 2 ,x' 3 co-moving with the rigid rotating object, 
where non-diagonal elements are null (due to the reflection symmetry of the mass 
distribution). Using the coordinate transformations from one system to another, one 
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can express the quadrupole tensor components in the x\,x 2 , £3 coordinate system as 
a function of the time-invariant D' a @ ones in the x[ , x' 2 , x' 3 system: 

Dn = D' n cos 2 (ujt) + D' 22 sin 2 (wt) 

D l2 = D 21 = X - sin(2wt)(D 11 - D' 22 ) (A.7) 

D22 = D' u sin 2 (wi) + D' 22 cos 2 (ojt) 
D 33 = D' 33 

From these expressions it is evident that D\\ = —D 22 and -D33 = 0, so that 
substituting in equations (|A.5|) and (|A.6[) one gets: 

hx = -^-(D' 11 -D 22 )sin(2cjt)cose (A.8) 

&+ = \^Zj( D 'ii - D 22) cos(2wi) (1 + cos 2 O ) (A.9) 

Considering that we expect to be observing the GRB on-axis (9q = in (|A.8|) 
and (|A.9[0 . the signal is circularly polarized. Equations (|A.8|) and (|A.9[) do apply also 
to the case of a binary system or to a bar-like structure, which are all thought to play 
a role in GRB progenitors [77] , 

An equivalent but useful way to expand the waveforms is in terms of I = 2 
pure-spin tensor harmonics. For a transverse, traceless tensor, in the quadrupole 
approximation, the only components that can enter are the basis states usually labeled 
by T E2 ' lm [86]. The expansion of the GW (transverse-traceless component of the 
metric perturbation) in this basis can be written as: 

hl T = -%A 2m T E2 - 2 ™ (A.10) 

where there is an implicit summation over m. Using the explicit representation of 
these harmonics in the orthonormal spherical coordinates (see [SB] for details), it is 
possible to show that: (i) in the case of a mass distribution having a quadrupole mass 
tensor of the form (j A.4|) . which corresponds to a system with an x-y plane reflection 
symmetry (i.e. D\ 3 — D 23 = 0), one has A 2 ±\ = 0; (ii) for a rigidly rotating ellipsoid, 
since equations (|A.7|) are valid (i.e. D\\ = —D22 and -D33 = 0), one has ^2±o = 0. 
Thus, we can say that GW emission from binary systems or rotating rigid ellipsoids 
are dominated by the I = |m| = 2 mode, for which the wave amplitude is maximized 
along the rotational axis. 



Appendix B. Energy radiated in GWs 

In what follows we give details on the procedure followed to determine the energy 
upper-limit for the sine-Gaussian waveform with Q = 5 at ~ 200 Hz. 
The energy radiated in GWs is computed as: 

r 3 d 2 f /"+ 00 • • dt 

Egw = t^gJ dn j x {h + {t) + hl{t)) rr z (R1) 
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where the integration over the solid angle should be performed while considering the 
emission pattern. If the signal power at the detectors is dominated by a frequency / , 
as is the case for the sine-Gaussian waveforms, the above formula is approximated as: 



c3d h(^ 2 fZ) I '<n [ +c °(hi(t) + hi(t)) 



16wG 

On the basis of Eqs. (|A.8|) - (|A.9| ), we have: 

3j2n 2f2 r i 



l + z 



(B.2) 



dt x 



-(1 + cos 2 Ofh 2 + Q {t) + cos 2 6h\ Q {t) 



(B.3) 
(B.4) 

where h+ t o(t) and /i x ,o(0 are given by the plus and cross components in Eq. (p~0|) - (fTTj) . 
Taking into account that h\ dt = h 2 + dt, we can write: 



Egw 



3^2 f2 i-l 

' d(cos6») 



4G(l + z) 
dt h\ Q (t) 



-(l + cos 2 ey + cos 



(B.5) 

The time integral is equal to (h^ s ) 2 /2 7 where h^ s has the values quoted in Table [T] 
Thus we write: 



Egw — (Kss) 



SG\2 c3 ^|27t 2 /q 



5G(l + z) 
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